Density Matrix Renormalization Group Study of the 
Disorder Line in the Quantum ANNNI Model 
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We apply Density Matrix Renormalization Group methods to study the phase diagram of the 
quantum ANNNI model in the region of low frustration where the ferromagnetic coupling is larger 
than the next-nearest-neighbor antiferromagnetic one. By Finite Size Scaling on lattices with up 
to 80 sites we locate precisely the transition line from the ferromagnetic phase to a paramagnetic 
phase without spatial modulation. We then measure and analyze the spin-spin correlation function 
in order to determine the disorder transition line where a modulation appears. We give strong 
numerical support to the conjecture that the Peschel- Emery one-dimensional line actually coincides 
with the disorder line. We also show that the critical exponent governing the vanishing of the 
modulation parameter at the disorder transition is /3 q = 1/2. 
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The ANNNI model is an axial Ising model with com- 
peting next-nearest-neighbor antiferromagnetic coupling 
in one direction. It is a paradigm for the study of compe- 
tition between magnetic ordering, frustration and ther- 
mal disordering effects. Its phase diagram displays in- 
deed a rich variety of phases. In the most realistic three- 
dimensional case, it describes several physical systems 
from magnetic materials like CeSb to binary alloys or 
dielectrics like NaNC>2 Q- In the more academic one- 
dimensional case, it is exactly solvable and several gen- 
eral properties can be rigorously proved about its phase 
diagram 0|. 

The two-dimensional case is nontrivial and not solv- 
able, but its phase structure is much simpler than the 3d 
case. The model is believed to display 5 phases 0: Fer- 
romagnetic TTTTj antiphase TTII, paramagnetic with or 
without modulation and floating phase with algebraically 
decaying spin correlations. This picture is supported 
by various analytical 0] and numerical [j| investigations 
based on a variety of approximations. However, lack- 
ing an exact solution, the precise location of the various 
transitions is not known beyond approximate treatments. 
Actually there is no rigorous proof of the existance of all 
the above phases. In particular the very existance of the 
floating phase has been recently under debate 0. 

To further simplify the analysis, the 2d case can 
be studied in the Hamiltonian limit which is a one- 
dimensional quantum spin S = 1/2 chain with next- 
nearest-neighbor coupling. The chain interacts with an 
external field playing the role of the temperature and 
triggering phase transitions. The Hamiltonian limit, also 
called TAM model (Transverse ANNNI) is very interest- 
ing in itself being a simple example where several com- 



plicated quantum phase transitions do occur with drastic 
changes in the qualitative features of the ground state • 

The accurate numerical study of the TAM model is 
challenging notwithstanding its relative simplicity. In 
this Letter we address an open conjecture concerning its 
disorder line by employing Density Matrix Renormaliza- 
tion Group methods. 

To illustrate the problem, we introdude the TAM 
Hamiltonian with open boundary conditions which reads 
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We shall present our results in terms of the adimcnsional 
ratios K = —J2/J1 and B / J\ which are the only param- 
eters that describe the properties of the ground state. 

The qualitative phase diagram of the TAM model is 
quite different in the two regions k < 1/2 or k > 1/2. 
For k < 1/2, the model is in a ferromagnetic phase at 
low magnetic field B. At B c \, a transition in the Ising 
class makes the ground state paramagnetic with exponen- 
tially decaying spin-spin correlation functions. Increasing 
further the external field we expect a new transition for 
B > B C 2 > B Cj i where the model is still gapped but with 
a correlation function whose exponential decay has also a 
spatial modulation. In this phase the asymptotic corre- 
lation function in the bulk is conveniently parametrized 
for large spin separation d by the functional form 

C zz {d) = {of af +d ) ~ c e- d / r cosfrqd + <p), (2) 

with r and q being functions of B and k. 
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The modulation parameter q(B, re) vanishes at B = 
-Be, 2 («0 with a certain exponent (3 q 

q(B,K)~A(B-B Ct2 (K)) l3 '>, as B \ B c , 2 (k). (3) 

The critical line B = -B Ci 2(re) is known as a disorder 
line 0, (see also [jj for a different definition) . 

The region re > 1/2 is much more complicated. At low 
B the ground state is in a so-called antiphase with typi- 
cal spin configuration ft II ' ' ' • Increasing the magnetic 
field one expects to observe a first transition to a disor- 
dered phase with algebraically decaying C zz (the floating 
phase) followed by a final transition to the asymptotic 
paramagnetic phase, i.e., the unique high temperature 
phase in the original 2d statistical model. The numerical 
data in this region are controversial and the size of the 
floating phase is not clear being possibly zero |6(. 

In this Letter, we shall be concerned with the re < 1/2 
region only. For simplicity, we shall denote this region by 
LFR (low frustration region). In the LFR, there is gen- 
eral consensus about the phase diagram, although only 
at the qualitative level, i.e., with large variations due to 
the various approximation employed in its calculation. 

Remarkably, the TAM model can be solved exactly 
on a critical line in the LFR called the Peschel-Emery 
one-dimensional line (ODL) [l(l |. The spin correlation 
decays exponentially on the ODL which is immersed in 
the paramagnetic phase. Little is known analytically off 
the ODL due the the very tricky nature of the solution. 
It is still a conjecture that the ODL is indeed the disorder 
line and that therefore 

B c>2 /J x =B^ E /J l = K-^-. (4) 

4re 

The conjecture is compatible with the numerical sim- 
ulations of the TAM model. However, the agreement 
B c ,2 — B PE is valid at not more than about 20 % accu- 
racy along the line. 

The aim of this Letter is precisely to give a numer- 
ical proof with good accuracy of this conjecture. As a 
byproduct we also determine the unknown exponent (3 q . 

A detailed analysis of the quantum ANNNI model can 
be found in ^lj. The accuracy of the results is poor 
because of the small considered lattices with less than 10 
sites. Another interesting approach is described in 
where an effective Hamiltonian is proposed allowing a 
considerable reduction of the Hilbert space. Systems up 
to 32 sites long can be treated, but the approximation is 
valid only near re = 1/2. 

A more recent numerical analysis of the LFR is 0] 
where the ferromagnetic-paramagnetic Ising transition is 
analyzed by combining Finite Size Scaling with exact di- 
agonalization of short chains with no more than 10 sites. 

Here we present a study of the model with higher ac- 
curacy and much larger lattices by means of the Density 
Matrix Renormalization Group (DMRG) algorithm [14j . 
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FIG. 1: Finite Size Scaling analysis of LAi(k, B) at re = 0.3. 



Nowadays, this method appears to be the natural choice 
for one-dimensional quantum spin chains. 

We have implemented the finite size version of the 
DMRG algorithm computing the two lowest levels i?o,i 
and the energy gap A = E\ — Eq . The algorithm results 
are very stable when more than 80 states are kept in the 
block Hamiltonians. In practice the numerical error on 
A is at the level of the machine precision. 

For several lattice sizes L of order 10 2 and various frus- 
tration ratios re, we have computed the scaled energy gap 
LA l (k, B). The crossing of the associated curves as a 
function of B at fixed re is a finite size estimate of the 
ferromagnetic critical field B^[ (re). As an example, we 
show in Fig. (Q) our results at re = 0.3. We have deter- 
mined the crossing point between the curves associated 
to a certain L and L + 10. We expect B^i - * 5 C) i(re) 
as L — ► co with algebraic corrections in 1/L [15|. We 
show in Fig. the finite size crossing field plotted as a 
function of x = 1/(L + 10) which is the most convenient 
variable to extrapolate our (L, L + 10) crossings. Indeed, 
the fitting function a + bx 2 + cx 3 gives a very good x 2 
of about 10~ n . The results for all the considered re are 
collected in Tab. {1} where we also show (when available) 
the analogous results from These are obtained on 

small lattices crossing L with a fixed L = 4. This is at 
most an estimate of B c l . Tab. |IJ reports also B c i which 
is obtained from the vanishing of the gap at second order 
in B |l(i|. and is defined by 

1 + 2k= J^ + k J^±). (5) 
Ji 2(1 + re) ^ J x ) y > 

In principle, it is possible to determine -B c ,i directly in 
the infinite size limit by using the infinite lattice version 
of the DMRG algorithm. We show in Fig. © the result 
of such a procedure at re = 0.4. The result for B c \ is 
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FIG. 2: Extrapolation of the Finite Size Scaling crossing es- 
timator . k = 0.3 and L = 10, 20, . . . , 80. 



fully consistent with the FSS analysis. Also, the expo- 
nent v = 1 which governs the vanishing of the mass gap 
is very clear at L = oo. For the other values of k we 
have preferred to avoid the infinite size algorithm since 
it is known that it can fail when the phase structure is 
complicated 0. From Tab. ® we see that the DMRG 
estimate B,P^ RG and the results of are globally simi- 
lar and slightly below the approximation -B Cj i, especially 
at large \k — 1/2 1. The value from [13J at k = 0.4 is some- 
what away from the common values of _B ( Pj VIRG and -B c ,i- 

After the determination of the ferromagnetic- 
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B c ,i/Ji [13] B cA /Ji 


-De, 2 /Jl 


B c , 2 /Ji [10j 


0.15 


0.73405(4) 


0.7327(2) 


0.74956 


1.5168(2) 


1.51667 


0.20 


0.6393(1) 


0.6407(4) 


0.65336 


1.0500(1) 


1.05 


0.25 


0.5403(3) 


0.5388(4) 


0.55051 


0.75001(2) 


0.75 


0.30 


0.43669(4) 


0.4368(2) 


0.44183 


0.53337(5) 


0.53333 


0.35 


0.32821(2) 


0.3298(3) 


0.32917 


0.36428(5) 


0.36429 


0.40 


0.216090(3) 0.2068(3) 


0.21548 


0.22498(2) 


0.225 



TABLE I: Comparison of the ferromagnetic and disorder crit- 
ical fields S Cj i_2/Ji. The DMRG columns are our data. _B c ,i 
is obtained from the vanishing of the gap at second order in 
the magnetic field. The data from |T3| have been interpolated 
at k = 0.15, 0.35. 



paramagnetic Ising transition, we studied B c ^{n) and 
the critical behavior of the modulation q. We have 
computed by the DMRG algorithm the spin correlation 
C zz on lattices large compared to the correlation length 
r appearing in Eq. In practice, L — 40 is enough in 
all the considered cases. The critical behaviour of q(B) 
is shown in Fig. Q for k — 0.3. The vanishing of q 2 is 
linear in B — B c i. The modulation parameter vanishes 
with exponent (3 q — 1/2. The critical field B C 2 coincides 
with the Peschel-Emery value with high accuracy. 

We have repeated the analysis for k — 0.15, 0.2, 
0.25, 0.35, 0.4, finding always a very good agreement. 
The agreement at small frustration is remarkable. From 
the point of view of the disorder line the next-nearest- 
neighbor coupling Ji is a singular perturbation with 
B c .2 —> oo in the isotropic Ising limit k —> 0. We remark 
that it is nontrivial to extend the calculation of off 
the one-dimensional line to proof rigorously that the one- 
dimensional line is the disorder line. Indeed, the only 
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FIG. 3: Finite Size Scaling analysis at k = 0.4. The left 
plot includes the data obtained with the infinite size DMRG 
algorithm and labeled L — oo. 



FIG. 4: B dependence of the squared modulation parameter 
q 2 at k — 0.3. The fit is performed on the leftmost points 
near the critical point. 
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[Eq. (5)] 
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FIG. 5: Phase diagram in the LFR showing the agreement 
between the disorder line computed with DMRG and the ex- 
act one-dimensional line by Peschel- Emery. See Tab. |lj for 
numerical values. 



analytic insight in this direction is the analysis in [l^ 
where the ODL is proved to be the disorder line, but 
only mapping the initial 5=1/2 spin chain into a dual 
spin T — 1/2 chain and taking the T — > oo limit. 

Our results for B C 2 are also summarized in Tab. (QQ) 
together with the Peschel-Emery value. In Fig. (JjJJ we 
plot the final phase diagram as determined by our DMRG 
simulations. 

In conclusion, we have shown that a DMRG analysis 
of the quantum ANNNI model provides strong numerical 
support to the conjecture that the Peschel-Emery ODL 
is actually the disorder line. Also, the critical exponent 
governing the vanishing of the modulation at the disorder 
transition is (3 q = 1/2. 

A natural extension of this work concerns the DMRG 
study of the region k > 1/2, which requires much larger 
lattices to analyze the slow algebraic decay of the spin 
correlation functions in the would-be floating phase. 

We acknowledge conversations with G. F. De Angelis, 
W. Selke, P. Fendley, and V. Rittenberg. 
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